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Abstract. A criterion for matrix failure of laminated composite plies in transverse tension and in-plane 
shear is developed by examining the mechanics of transverse matrix crack growth. Matrix cracks are 
assumed to initiate from manufacturing defects and can propagate within planes parallel to the fiber 
direction and normal to the ply mid-plane. Fracture mechanics models of cracks in unidirectional laminates, 
embedded plies and outer plies are developed to determine the onset and direction of propagation for 
unstable crack growth. The models for each ply configuration relate ply thickness and ply toughness to the 
corresponding in-situ ply strength. Calculated results for several materials are shown to correlate well with 
experimental results. 


1. INRODU C TION 

Strength-based failure criteria are commonly 
used with the finite element method to predict 
failure events in composite structures. A large 
number of continuum-based criteria have been 
derived to relate internal stresses and 
experimental measures of material strength to the 
onset of failure. The most typical among the 
proposed criteria is Hashin’s criterion [7] which 
assumes a quadratic interaction between the 
stresses acting on the fracture plane. 

The recent World Wide Failure Exercise 
(WWFE) conceived and conducted by Hinton 
and Soden [2] provides a good assessment of the 
status of currently available theoretical methods 
for predicting material failure in fiber reinforced 
polymer composites. Comparison of the 
predictions by the WWFE participants with 
experimental results indicates that even when 
analyzing simple laminates that have been 
studied extensively over the past 40 years, the 
predictions of most theories differ significantly 
from the experimental observations. 

The uncertainty in the prediction of initiation and 
progression of damage in composites has led to 
the undertaking of an effort at the NASA 
Langley Research Center to revisit existing 
failure theories, assess their capabilities, and to 
develop new theories where necessary. As a 
result of that effort, a set of six nonempirical 
criteria for predicting failure of unidirectional 
FRP laminates were proposed [ 3 ]. The objective 


of this paper is to present a more detailed 
examination of the failure criterion for matrix 
cracking presented in Ref. [ 3 ]. 

2. MATRIX CRACKING AND IN-SITU 
STRENGTHS 

Transverse matrix cracking is often considered a 
benign mode of failure because it normally 
causes such a small reduction in the overall 
stiffness of a structure that it is difficult to detect 
during a test. However, transverse matrix cracks 
also provide the primary leakage path for gases 
in pressurized vessels. Leakage is a phenomenon 
that is receiving considerable interest after the 
cancellation of the NASA X-33 launch vehicle 
program [ 4 \, where it was found that liquid 
hydrogen had leaked at cryogenic temperatures 
and, once the temperature increased, caused 
major delaminations by cryopumping. 

To predict matrix cracking in a laminate 
subjected to in-plane shear and transverse tensile 
stresses, a failure criterion must be capable of 
calculating the “in-situ” strengths. The in-situ 
effect, originally detected in Parvizi’s [5] tensile 
tests of cross-ply glass fiber reinforced plastics, 
is characterized by higher transverse tensile and 
shear strengths of a ply when it is constrained by 
plies with different fiber orientations in a 
laminate, compared with the strength of the same 
ply in a unidirectional laminate. The in-situ 
strength also depends on the number of plies 
clustered together, and on the fiber orientation of 
the constraining plies. 
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The orientation of the constraining plies and the 
number of plies clustered together also affect the 
crack density and the stiffness reduction of the 
cracked ply. Wang’s \6\ tests of [0/90n/0] 
(n=l,2,3,4) carbon/epoxy laminates have shown 
higher crack densities for thinner 90° layers. The 
reduction of the elastic properties of a cracked 
ply is normally predicted using elastic analyses 
of cracked plies [7, 8], or Continuum Damage 
Models [9-12]. 

The in-situ effect is illustrated in Fig. 1, where 
the relation between the in-situ transverse tensile 
strength and the total thickness of the 90° plies 
clustered together is represented. 

Accurate in-situ strengths are necessary for any 
stress-based failure criterion for matrix cracking 
in constrained plies. Both experimental [6, 13, 
14] and analytical methods [7, 15, 16] have been 
proposed to determine the in-situ strengths. In 
the following sections, the in-situ strengths are 
calculated using fracture mechanics solutions for 
the propagation of cracks in constrained plies. 



Fig. 1 - Transverse tensile strength as a 
function of number of plies clustered 
together, with models from Dvorak [75] 
based on experimental data from Wang \6\. 


Fracture Mechanics Analysis of a Cracked Ply 

The failure criterion for predicting matrix 
cracking in a ply subjected to in-plane shear and 
transverse tension proposed here is based on the 
fracture mechanics analysis of a slit crack in a 


ply, as proposed by Dvorak and Laws. [75] A slit 
crack represents the macroscopic effect of 
distributed matrix-fiber debonds that occur at the 
micromechanical level [6], as shown in Fig. 2. 
Slit cracks are the result of manufacturing 
defects or residual thermal stresses caused by 
different coefficients of thermal contraction 
between the matrix and the fibers. Before the slit 
crack grows unstably, it is idealized as lying on 
the 1-3 plane, as represented in Fig. 3. The crack 
has a length 2a 0 across the thickness of a ply, t. 



Fig. 2 - Fiber-matrix debonds coalesce to 
form a longer transverse crack. Ref. [7 7] 


Fig. 3 - Slit crack geometry (after 
Dvorak[75]). 


Plane stress conditions are assumed. The 
transverse tensile stress c 2 2 is associated with 
mode I loading, whereas the in-plane shear stress 
r 12 is associated with mode II loading. The crack 
represented in Fig. 3 can grow in the 1 
(longitudinal, L) direction, in the 3 (transverse, 
T) direction, or in both directions. 

The components of the energy release rate for 
the crack geometry represented in Fig. 3 were 
determined by Dvorak and Laws [75]. For 
mixed-mode loading, the energy release rate for 
crack growth in the T and L directions, G(T) and 
G(L j, respectively, are given by: 
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G (T) = ^l (4 a ° 22 <4 + tjI a ° 44 4 ) 

( 1 ) 

G(L) = ZjS- (#/ A° 22 <t 22 + 4 A° 44 4 ) 


G,(T) = ^A‘„(r;7 

(4) 

G nc (T) = ^KMsJ 


where it can be observed that the energy release 
rate G(L) for longitudinal propagation is a 
function of the transverse slit size a 0 and that it 
is not a function of the slit length in the 
longitudinal direction a q . 

The parameters ij 1 , i=I t II in Eq. 1 are the stress 
intensity reduction coefficients for propagation 
in the transverse direction, and the parameters 
i=I t II are the reduction coefficients for 
propagation in the longitudinal direction. These 
coefficients account for the constraining effects 
of the adjoining layers on crack propagation: the 
coefficients are nearly equal to unity when 
2a 0 «t , and are less than unity when a^t. 
Experimental results [14\ have shown an 
increase in the in-situ transverse tensile strength 
of / 1 0 90 n / s . 0=0°, 30°, 60°, laminates for 

increasing stiffness of adjoining sublaminates 
± 0 : This implies that the value of the parameter 
rfi decreases with increasing stiffness of 
adjoining sublaminates. Considering that a 
transverse crack can promote delamination 
between the plies, Dvorak and Laws [75] 
suggested that the effective value of rj, can be 
larger than obtained from the analysis of cracks 
terminating at the interface, and suggested the 
use of rji = £ =7. 

The parameters A° /; are calculated as [75] 


A = 2 

22 


V^2 


v 2 A 
U 2\ 


A 0 , 


( 2 ) 


The mode I and mode II components of the 
energy release rate can be obtained for the T- 
direction using Eq. (12) with 77 * = 7 


G,(r) = ^A> 2 2 2 

G n (T) = Zf- A°„4 


( 3 ) 


The corresponding components of the fracture 
toughness are given as: 


where Y t T s and Sf s are the in-situ transverse 
tensile and shear strengths, respectively. 

For propagation in the longitudinal direction, the 
mode I and mode II components of the energy 
release rate are: 


n st \ — 71 a 0 a 0 ^2 
G 1 \L) — A 22 <t 22 

( 5 ) 

a: / T \ ^ ^0 A 0 ^*-2 

Gn(L) — ^ A 44 T u 

and the components of the fracture toughness 
are: 


( 6 ) 

G nc (L) = ^ KM)' 


Failure Criterion for Matrix Cracking 

Having obtained expressions for the components 
of the energy release rate and fracture toughness, 
a failure criterion can be applied to predict the 
propagation of the slit crack represented in Fig. 
3. Under the presence of both in-plane shear and 
transverse tension, the critical energy release rate 
G c depends on the combined effect of all 
microscopic energy absorbing mechanisms such 
as the creation of new fracture surface. Relying 
on microscopic examinations of the fracture 
surface, Hahn [78] observed that the fracture 
surface topography strongly depends on the type 
of loading. With increasing proportion of the 
stress intensity factor K Ih more hackles are 
observed in the matrix, thereby indicating more 
energy absorption associated with crack 
extension. Therefore, Hahn proposed a mixed 
mode criterion written as a first-order 
polynomial in the stress intensity factors Kj and 
K n . Written in terms of the mode I and mode II 
energy release rates, the Halm criterion is 


(i-*l 


Gj(f) j Gj( i) | G n (i) l 

G Ic (i) G Ic (i) G IIc (i) 


i = T,L 


( 7 ) 
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where the material constant g is defined 
identically from either Eq. 4 or 6 as: 



A°,m 


( y T ^ 2 


\S is j 


( 8 ) 


A failure index for matrix tension can be 
expressed in terms of the ply stresses and in-situ 
strengths Yj s and S t L s by substituting either Eqs. 
3-4 or 5-6 into the criterion in Eq. 7 to yield: 


However, such tests demand complex set-ups, 
especially for the determination of the in-situ 
shear strength, and are expensive to perform. It is 
therefore desirable to calculate the in-situ 
strengths from material properties measurable 
from standard test specimens. In the following 
sections, various ply configurations are 
examined and fracture mechanics idealizations 
are applied to relate the stresses at failure to the 
ply toughnesses. 




\ 2 


y 

J 




= 1 


( 9 ) 


The criterion presented in Eq. 9, with both linear 
and quadratic terms of the transverse normal 
stress and a quadratic term of the in-plane shear 
stress, is similar to the criteria proposed by Hahn 
[18], Liu [19] (for transverse tension and in- 
plane shear), and Puck [20]. In addition, if g=l, 
Eq. 7 reverts to the linear version of the criterion 
proposed by Wu and Reuter [21] for the 
propagation of delamination in laminated 
composites: 



( 10 ) 


Furthermore, using g=l, Eq. 9 reverts to the 
well-known Hashin criterion [1] for transverse 
matrix cracking under both in-plane shear and 
transverse tension, where the ply strengths are 
replaced by the in-situ strengths: 


FT, 


Y 1 

\Fs j 




<1 


( 11 ) 


Calculation of In-Situ Strengths 

The failure criterion proposed in Eq. 9 can only 
be as accurate as the values of the in-situ 
strengths Y? and that are used. These 
strengths represent the averaged ply stresses at 
which transverse cracks grow unstably. As was 
described earlier, these strength values depend 
on ply configuration parameters such as ply 
thickness and ply toughness. 

It is possible to perform experimental tests to 
measure the in-situ strengths of laminates. 


Failure of Thick Embedded Plies 

A thick ply is defined as one in which the length 
of the slit crack is much smaller than the ply 
thickness, 2 a 0 « t , as illustrated in Fig. 4. The 
minimum thickness for a thick ply depends on 
the material used. For E-glass/epoxy and 
carbon/epoxy laminates, Dvorak and Laws [15] 
calculated the transition thickness between a thin 
and a thick ply to be approximately 0.7 mm, or 
about 5 to 6 plies. 



‘2a n 


© 



Fig. 4 - Geometry of slit crack in a thick 
embedded ply subjected to tension and shear 
loads. 


For the geometry represented in Fig. 4, the crack 
can grow in the transverse or in the longitudinal 
direction. Comparing Eqs. 3 and 5, however, 
indicates that the energy release rate for the 
crack slit is twice as large in the transverse 
direction as it is in the longitudinal direction. 
Since Eq. 3 also indicates that the energy release 
rate is proportional to the crack length 2 a Q , the 
crack will grow unstably in the transverse 
direction. Once the crack reaches the 

constraining plies, it can propagate in the 
longitudinal direction, as well as induce a 
delamination. 

The condition for crack propagation under 
biaxial loads is calculated using Eq. 9, where the 
in-situ strengths are obtained from the 

expressions of the fracture toughness in Eqs. 4: 
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2GjJT) 

, — ( 12 ) 

s l = 2 G IIc (T) 

is i* a oK 4 

It can be observed from Eqs. 12 that the in-situ 
strengths of thick plies Y? and are functions 
of the toughnesses G lc (T) and G IIc (T ) of the 
material and the size of the material flaw, 2 a 0 . 
Therefore, the in-situ strengths for thick plies are 
independent of the ply thickness, as has been 
observed by Dvorak [75] andLeguillon [22], and 
as was shown in Fig. 1 . 

Failure of Thin Embedded Plies 

Thin plies are defined as having a thickness 
smaller than the typical defect, t < 2a 0 , so the 
slit crack extends across the entire thickness t of 
the ply, as represented in Fig. 5. 



Fig. 5 - Geometry of slit crack in a thin 
embedded ply. 


experimentally obtained in-situ strengths of both 
thick and thin 90° plies in [0/90/0] laminates, as 
was shown in Fig. 1. 

Failure of Unidirectional Laminates 

The fracture of a unidirectional specimen is 
taken as a particular case of a thick ply [75]. The 
defect size is 2 a 0 , as for thick embedded plies. 
However, in the absence of constraining plies, 
the critical initial slit crack is located at the 
surface of the laminate. For tensile loading, the 
crack can be located at the edge of the laminate, 
which increases the energy release rate when 
compared with a central crack. In the case of 
shear loading, there is no free edge, so the crack 
is a central crack, as shown in Fig. 6. 



Crack propagation for unidirectional specimens 
subjected to tension or shear loads is obtained 
from the classic solution of the free edge crack 
[25]. 


In the case of thin plies, crack defects can only 
grow in the longitudinal (L) direction, or trigger 
a delamination between the plies. The in-situ 
strengths for thin plies can be calculated from the 
components of the fracture toughness as: 


Thin: 


*G & (L) 
K t A 22 


g GnSPj 

71 1 A° 44 


(13) 


where it can be observed that the in-situ 
strengths are inversely proportional to yft . The 
toughnesses G lc (L) and G Uc (L ) can be 
assumed to be the values measured by standard 
Fracture Mechanics tests, such as the double 
cantilever beam test (DCB) for mode I and the 
end notch flexure test (ENF) for mode II. Using 
Eqs. 12 and 13, Dvorak and Laws [75] obtained 
a good correlation between the predicted and 


Gi c (T) = \.n 1 Ka 0 \ 0 21 (Y T Y 

G Ilc (T) = xa li A° u (s L ) 2 

where Y T and S L are the material tensile and 
shear strengths as measured from unidirectional 
laminate tests. Note that the in-situ strengths for 
thick plies can be obtained as a function of the 
strength of unidirectional laminates by 
substituting Eqs. 14 into 12: 


Y* = l.l2y[2Y T 
Thick: _ 

=V2 s L 


(15) 


The failure criterion for unidirectional plies 
under in-plane shear and transverse tension is 
represented in Eq. 9. The toughness ratio g for a 
thick laminate can also be calculated in terms of 
the unidirectional properties by substituting Eqs. 
14 into Eq. 8 to yield: 
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g = 



1.12 2 


A 0 ,. 


/ 


V 



(16) 


crack growth in the longitudinal direction, for 
both mode I and mode II loading, is also required 
for the determination of the in-situ strengths. 


Failure of Outer Plies 

Outer plies are taken as a special case of thin 
plies, as represented in Fig. 7. Following the 
procedure presented for a unidirectional 
laminate, it is considered that the stress intensity 
factor of an outer ply is larger than the stress 
intensity factor of a thin embedded ply. The 
relation between the stress intensity factors given 
by classical Fracture Mechanics solutions for 
free-edge cracks. 




Fig. 7 - Geometry of slit crack in a thin 
outer ply. 


In the case of thin outer plies, crack defects can 
only grow in the longitudinal (L) direction, or 
trigger a delamination between the plies. The in- 
situ strengths of thin outer plies can be calculated 
from the components of the fracture toughness 
as: 


The in-situ strength factors for thick embedded 
plies are given in Eqs. 15. The in-situ strength 
factors for thin embedded plies and for thin outer 
plies are calculated by substituting Eqs. 13 and 
17 into Eq. 18. The in-situ factors for the 
different configurations are shown in Table 1. 


Table 1 - In-situ Factors in Modes I and II 



Mode I 

Mode II 

Thick 

embedded 

ply 

(Oj= I. 12 V 2 

CO n = 

F 

Thin 

embedded 

ply 

F 

G IC (L ) 

CO -F I 

\G IIC (L) 

7 r T ] 

7rA° 22 t 

77 S l ] 

M.V 

Outer ply 

_ 1.78 

G IC (L) 

2 

/» — 

IGjjc(L) 

U)j 

Y r 1 

7rA ° 22 t 

77 ] 

1 y 


The in-situ factors for the transverse tensile 
strength of T300/934 CFRP are presented in 
Fig. 8. The material properties were obtained by 
Dvorak and Laws [75]. 


Outer: 


= 1.79 r Ic f 



\ n tA 22 


«i 

(17) 

4- 

=2 



i ntK\, 


3.5- 


In-Situ Strength Factors 

The in-situ matrix tensile and shear strength 
properties Y* and Sf s for composite plies are 
often expressed in terms of the measured 
unidirectional strengths Y T and S L using in- 
situ strength factors C0j and 0) n that reflect the 
ply configuration. In general, the in-situ 
strengths are 


T300/934 CFRP 

\ G lc (L)=220J/m 2 

\ Y T =40.7MPa 


\Thin embedded ply 


2 - 


1.5- 


,Thin outer ply v 


Thick ply 


St=co n S L 


(18) 


UD laminate 

0.2 0.3 0.4 0.5 0.6 07 0.8 0.9 

Ply thickness (mm) 


Fig. 8 - Predicted Mode I in-situ factors for 
T300/934 CFRP. 


The unidirectional transverse tensile and in-plane 
shear strengths are determined from standard 
experimental tests. The fracture toughness for 
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3. COMPARISON WITH EXPERIMENTS 

A comparison of the strength of an AS4-55A 
carbon fiber reinforced ply predicted using Eq. 9 
and published experimental results by Swanson 
[24] is shown in Fig. 9. The calculations were 
based on the material properties shown in Table 
2. It can be observed that the Hasliin criterion 
overpredicts the strength when shear is the 
dominant load. The results predicted using Eq. 9 
are in overall good agreement with the 
experimental results. 



Fig. 9 - Comparison between calculated and 
experimental data for AS4/55A 
carbon/epoxy. (Swanson, Ref. [24]) 


Table 2 - Properties of AS-4 55A in MPa. [24] 

E„ E 22 G i2 7 V 

126,000 11 ,000 6,600 0.28 26.7 51.8 


The strength of a glass/epoxy Scotch-Ply Type 
1002 ply predicted using Eq. 9 is compared in 
Fig. 10 with published experimental results by 
Voloshin [25]. The calculations were based on 
the material properties shown in Table 3. It can 
be observed that the Hasliin criterion 
overpredicts the strength when shear is the 
dominant load. The results predicted using Eq. 9 
are in good agreement with the experimental 
results. 


Table 3 

- Properties of Scotch-ply E-glass 
MPa. [25] 

in 


E n 

E 2 2 Cj 2 

Y T 

V 

53,500 

17,700 5,830 

0.28 20.6 

37.5 



Fig. 10 - Comparison between calculated 
and experimental data for glass/epoxy 
Scotch-Ply Type 1002. (Voloshin [25]) 


4. DISCUSSION 

Accurate failure criteria are necessary to assess 
the structural integrity of laminated composite 
structures. Transverse matrix cracking is a failure 
mode that is often considered benign and not 
given detailed consideration. However, 
transverse matrix cracks are the primary leakage 
path for gases in pressurized tanks. In the present 
paper, matrix cracking under combined 
transverse tension and inplane shear is discussed. 
In addition, the issue of the in-situ strength is 
examined by developing fracture mechanics 
idealizations of different ply configurations. The 
proposed fracture mechanics models relate the 
effective ply strength to the ply toughnesses, 
while accounting for the constraining effects of 
adjacent plies. 

Failure envelopes were constructed for a 
graphite/epoxy material and a glass/epoxy 
material, and the predicted results indicate good 
correlation with experimental results. 
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